Non-equilibrium transport through a vertical quantum dot 
in the absence of spin-flip energy relaxation 
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We investigate non-equilibrium transport in the absence of spin-flip energy relaxation in a few- 
electron quantum dot artificial atom. Novel non-equilibrium tunneling processes involving high- 
spin states which cannot be excited from the ground state because of spin-blockade, and other 
processes involving more than two charge states are observed. These processes cannot be explained 
by orthodox Coulomb blockade theory. The absence of effective spin relaxation induces considerable 
fluctuation of the spin, charge, and total energy of the quantum dot. Although these features are 
revealed clearly by pulse excitation measurements, they are also observed in conventional dc current 
characteristics of quantum dots. 



A quantum dot (QD) is a small conducting island in 
which electrons occupy discrete energy states [Q. En- 
ergy relaxation from an excited state (ES) to a ground 
state (GS) inside a QD is significantly suppressed if a 
spin flip is required [||, ^, ^ . Our measurements on QDs 
in the Coulomb blockade (CB) regime indicate an ex- 
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long spin- flip energy relaxation time, Tg- 



> 1 



which is much longer than the interval of tun- 
neling events (V~^ = 1 - 100 ns for a typical tunneling 
current of 1 - 100 pA, where V is the tunneling rate) 
as well as the momentum relaxation time, Tmo ~ 1 ns 
dl, 0, H]. In the absence of efficient energy relaxation, ex- 
cess energy remains in the system. This opens up novel 
non-equilibrium transport channels, which we describe in 
this Letter. 

We start from the orthodox theory which describes CB 
and single electron tunneling (SET). The total energy, 
U{N)^ of the system, in which an island containing N 
electrons is affected by a gate voltage, Vg, via a capaci- 
tance, Cn, is 
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[|l|, ^. The first term is the electrostatic energy approx- 
imated by a constant Coulomb interaction. Inside the 
parenthesis is the sum of the charge on the dot, the in- 
duced charge by the gate, and an offset charge, qo- Cy. 
is the total capacitance of the dot. The second term, 
Eint{N)^ is the sum of the energies of the occupied N 
electron levels, measured relative to the Fermi energy of 
the leads, accounting for the internal degrees of freedom 
of the QD. Other corrections to many-body interactions 
are also included in Eint- In the orthodox theory, orig- 
inally considered for a continuum density of dot states, 
the second term is neglected because the QD is assumed 



to relax quickly to the minimum energy, El 
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which is 



independent of N. In this SET scheme, an electron that 
has entered the island leaves before another electron is 



allowed to enter when U{No) = U^Nq + 1) as shown in 
Fig. 1(a). This situation is maintained unless the exci- 
tation energy exceeds the charging energy, Ec = e^/C^]. 
In this picture spin is neglected and A/'-electron GSs are 
only considered. 
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FIG. 1: (a) Total energy U{N) of a small classical island 
when SET allows N to fluctuate between A^o and A^o + 1. (b) 
U{N,S) of QD-1 at B = 3 T and Vg - -1.7 V. The arrows 
indicate allowed tunneling transitions with excitation energy 
< Ea (see text), (c) Schematic setup for the pulse measure- 
ment on a vertical QD. The circular pillar is of diameter 0.54 
/im (0.50 /im) for QD-1 (QD-2). See Ref. @ for details. All 
the measurements are performed at a temperature ^ 0.1 K. 
The magnetic field is applied parallel to the current, (d) The 
time-dependent gate voltage Voit). 

For a QD in which energy quantization and many-body 
interactions are significant, we must consider the discrete 
energy of the dot, Eint{N^ S^M)^ which is characterized 
by total spin, and total angular momentum, M [Q. 
We focus on the regime Tspin ^ F"^ ^ Tmo-, where spin- 
flip energy relaxation is effectively absent. This is the 



2 



typical condition for a dot weakly coupled to the leads 
- the coupling to the leads is still strong enough to give 
a measurable current (eP > 1 fA), but weak enough to 
prevent cotunneling current. // the QD is excited to any 
N- electron state with a different total spin from that of 
the N- electron GS, the ES cannot always relax to the GS 
before the QD undergoes a tunneling transition to an- 
other ± 1 electron state. Successive tunneling tran- 
sitions force the QD into highly nonequilibrium config- 
urations. Nonetheless, there is always a selection rule 
for the tunneling transition to be satisfied: each tun- 
neling transition changes N by one and S by one half. 
Tunneling transitions which change S by more than one 
half should be blocked (spin blockade) Figure 
1(b) shows a particular U{N^S) diagram, which can ac- 
tually be realized in our QD (see below). Long lived ESs 
are now included, and the different spin states have dif- 
ferent energies because of direct Coulomb and exchange 
interactions The allowed tunneling transi- 

tions indicated by the arrows require an excitation energy 
smaller than the addition energy, Ea. All these transi- 
tions can cause the dot state (A", and U) to fiuctuate 
dramatically. 

In order to investigate highly non-equilibrium trans- 
port, we employ a pulse excitation technique Q, which 
generates only transient current associated with long- 
lived spin states^ on two vertical QDs (QD-1 and QD- 
2) |l3|, Electrons are confined laterally by an 

approximate two-dimensional harmonic potential (con- 
finement energy hcjo ^ 4 meV for QD-1 and ~ 2.5 meV 
for QD-2). These samples show qualitatively the same 
behavior. The A^-dependent addition energy, Ea = 2 - 5 
meV, of QD-1 clearly reveals a shell structure |l^. A", 
S and M can be identified from the magnetic field, 5, 
dependence of the SET current spectrum. Zeeman split- 
ting is not resolved so we neglected it. A square pulse 
of amplitude Vp, combined with the static gate voltage, 
Vg, is applied to the gate electrode (g) shown in Fig. 1 
(c) 0. The time dependent gate voltage, Vcit), illus- 
trated in Fig. 1(d), is Vg = Vh = ^ ^Yp during the 
high phase of the pulse, and Vg = Vi = Vg — ^Vp during 
the low phase. A small dc bias voltage Vb = 0.15 mV 
is applied so that electrons are injected from the bottom 
contact b (F^^ ~ 10 ns), and escape to the top contact t 
(F^^ - 100 ns). The averaged dc current, /, is measured 
during the Fo-pulse irradiation. 

Figures 2(a) and (b) show pulse excited current spectra 
of QD-1 for A" = - 5 at 5 = 3 T. A current peak 
initially observed in the absence of the Fc-pulse is always 
split into two peaks of equal height, when the ^G-pulse is 
applied. Weaker additional peaks indicated by the arrows 
are due to transient current through ESs [||. The pulse 
length {th = ti = 300 ns) is chosen to be much longer than 
Tmo, but much shorter than Tspin- In this case, transient 
current through an A^-electron ES is only observed when 
the ES has a total spin different from that of any other 
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FIG. 2: (color), (a) Pulse excited current / vs. gate voltage, 
Vg, for QD-1 measured at B = 3 T and Vp ^ SO mV [dashed 
line in (b)]. Each horizontal bar indicates the Vg region where 
N can fluctuate, (b) Pulse excited current plot at B = 3 T 
built up from / vs. Vg. The color represents current ampli- 
tude [white (0 pA) - red (0.25 pA) - black (> 0.5 pA)]. N is 
fixed in each triangular region along the bottom due to CB. 
(c) Representation of current peaks (solid lines) observed in 
(b). The transient current peaks are marked by a, a, b, b^, 
etc. The prime (non-prime) indicates transport during the 
low (high) phase of the pulse. These peaks are extrapolated 
to = V as shown by dashed lines. The tunneling transi- 
tions between (N — 1, S ^ ^) and (A, S) states are also given. 
The red lines indicate current peaks attributed to the novel 
DET process, while cross-hatched regions indicate where nor- 
mal DET occurs, (d) Pulse excited current plot [white (0 pA) 
- red (0.15 pA) - black (> 0.3 pA)] between A = 2 and 3 at 
B — 3.7 T. The inset schematically shows the termination of 
line c by line b^ (e) The / vs. Vg data at various Vp values. 
The curves are offset for clarity. 

Examples of transient current excited by the VG-puhe 
for A" = 1, 2, and 3, are depicted in Fig. 3(a) (c) and (e), 
where U (A', S) is plotted against Vg- The normal GS-GS 
tunneling current peaks occur at level crossings labeled 
o, whilst transient current appearing at higher energy 
level crossings are denoted by different symbols. Figure 
3(a) shows U for N = 1 and 2, where a GS-ES tunneling 
transition labeled • takes place. Consider that a Vc-pulse 
is applied as indicated. Suppose the system is in the GS 
(A", S) = (1, ^) during the low-phase of the pulse (Vg = 
Vi). When the high-phase of the pulse is applied {Vg = 
14), the system changes non-adiabatically from (1, ^) at 
Vg = Vi to {1, ^) dit Vg = Vh (indicated by arrow i), since 
the rise time of the Fc-pulse (1-2 ns) is shorter than 
r-\ {Ttot = F, + F5). If Vh is tuned so that [7(1, ^) = 
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U (2, 1), transient current can flow at the crossing marked 
•. Only a few electrons can tunnel through the ES(2, 1), 
before the inelastic tunneling transition changes the state 
from (1, ^) to (2,0) This transient current appears if 
the pulse excitation energy, aVp, exceeds the level spacing 
A2 = /7(2, 1) - /7(2, 0) as indicated, where a = d[U{N) - 
U {N — l)]/dVG is the conversion factor from gate voltage 
to energy. The condition for the transient current to flow 
is given by the thick line marked • in the aVp - Vg plane 
of Fig. 3(b). This current is clearly seen as peak a in 
Fig. 2(b) 




FIG. 3: (a), (c), and (e): The gate voltage, Vcit), depen- 
dence of the total energy. The electron number N and total 
spin S are indicated by {N, S) . The arrows indicate the ex- 
citation/relaxation processes discussed in the text. The sym- 
bols indicate transport processes (o for GS-GS tunneling, • 
for GS-ES tunneling, for spin-blockade, for ES-ES tun- 
neling, and -k for novel DET). (b), (d), and (f): Correspond- 
ing conditions for stable GS-GS tunneling current (thin lines) 
and transient current (thick lines). The horizontal axis is the 
static gate voltage, Vg, while the vertical axis is the excitation 
energy of the Vc-pulse (= aVp). 

Next, we consider striking features for which the con- 
ventional SET scheme collapses. Consider the U {N^ S) 
for the = 1, 2, and 3 states in Fig. 3(c). Suppose a 
Vb-pulse, aVp < £^0, is applied as shown. If the state is 
initially (1, |) at Vq = VJ, /7(1, ^) is raised to a higher 
energy at Vg = Vh (arrow i), and subsequent inelastic 
tunneling results in either the (2,0) or (2,1) state. If it is 
(2, 1) (arrow ii), it is now possible for another electron to 
be transported at the crossing [7(2, 1) = [7(3, \) marked 
by -k. In this case transient transport persists until the 



inelastic tunneling transition from (3, \) to (2,0) occurs, 
because spin-flip energy relaxation from (2,1) to (2,0) is 
absent. The QD returns back to the initial state (1, ^) 
when the pulse is switched to Vi . This new tunneling pro- 
cess involves three charge states, and thus we consider it 
to be novel double electron tunneling (DET). Conven- 
tional DET within the orthodox CB theory appears only 
for aVp > Ea [|l|, |9|, while the novel DET takes place 
even for aVp < Ea. This process should appear on the 
thick line b in the aVp - Vg plane of Fig. 3(d). 

The novel DET is clearly seen in the experimental data 
of Fig. 2(b) as marked by b. Note that this peak cannot 
be due to regular SET involving higher lying N = 2 ESs. 
The spin singlet (TV = 2 GS at 5 = 3 T) and the spin 
triplet [peak a in Fig. 2(b)] are the only possible spin 
conflgurations. So no other ESs with different total spin 
can appear in the pulse measurement. Therefore, the 
pulse measurement successfully allows us to identify the 
extra peak b as a novel DET feature. Note that we can 
also conflrm that peak b is due to novel DET from the 
^-dependence of the peak position ]T6[ . 

Transient current can also appear during the low phase 
of the pulse. If a pulse is applied to excite the QD along 
(3, \ ) in the direction of decreasing Vg [dashed arrow 
in Fig. 3(c)], GS-ES tunneling between (3, \) and (2, 1) 
should appear along line b' in Fig. 3(d). Similarly, novel 
DET appears at the crossing /7(2,1) = /7(1,^) [corre- 
sponding to feature a' in Fig. 3(d)]. Note that the ex- 
trapolated lines a and dJ (b and b') meet at zero excita- 
tion energy. These features are clearly seen in Fig. 2(b). 

We now discuss spin-blockade and associated non- 
equilibrium transport. Figure 3(e) is the energy diagram 
for = 2 and 3 - the flrst case where spin-blockade 
appears. The direct transition between the GS (2,0) 
and ES (3, f ) is spin-blockaded (marked (g)). However, 
non-equilibrium tunneling between the ES (2, 1) and ES 
(3, |) (marked 0) is allowed if a pulse is applied as in- 
dicated. Even though the QD is initially in the GS 
(2,0) at Vg = Vj, it flrst changes to (3, \) at Vg = Vh 
(arrows i and ii), then to either (2,1) or (2,0) during 
the next low-phase [arrows iii and iv for (2,1)]. If it 
is (2,1), ES-ES tunneling occurs at /7(3, f ) = /7(2, 1) 
when the high-phase is restored (arrow v). This kind 
of tunneling mechanism leads to the fluctuation of the 
total spin from to |. Note again that this non- 
equilibrium fluctuation is attributed to the absence of 
effective spin relaxation. This ES-ES tunneling process 
requires complex excitation. The excitation energy must 
be greater than the corresponding two level spacings: 
aVp > A2 = U(2,l) - [7(2, 0) for the excitation indi- 
cated by arrow iii, and aVp > A3 = /7(3, |) — /7(3, |) 
for the excitation indicated by arrow v. The condition 
for ES-ES tunneling is marked by line c (0) in Fig. 3(f). 
This line is terminated by line b' (•) for GS-ES tunneling 
at Vg = VJ, so line c does not reach the line o for GS-GS 
tunneling. This termination is the signature of an ES-ES 
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tunneling process which does not involve any GS. 

We can identify these features in our QD. Peak c in Fig. 
2(d) and (e) now measured at 5 = 3.7 T is attributed to 
ES(2, 1) - ES(3, |) tunneling from the 5-field dependence 
of the peak position. This peak is clearly terminated by 
peak b', which is similarly assigned to ES(2, 1) - GS(3, |) 
tunneling. However, no measurable current (< 10 fA) 
is seen for spin-blockaded tunneling between (2,0) and 
(3, 1) in the region Vp = 50 - 100 mV and 5 = 3.5 - 4.1 
T (not shown) where spin blockade is expected. These 
observation are consistent with the above explanation. 
Note peak c is too weak to see at 5 = 3 T [Fig. 2(b)]. 
Nonetheless, the following are very clear in Fig. 2(b) 
for AT = 3 and 4: ES (3, f ) - ES (4,1) tunneling line 
d' is terminated by GS (3, ^) - ES (4, 1) tunneling line 
e, and the termination is associated with spin-blockaded 
ES(3, f) - GS(4,0) tunneling (no signal). 

Other complicated tunneling processes involving both 
novel DET and ES-ES tunneling are observed. Peak d 
(c^) in Fig. 2(b), which is terminated by peak b' (e), 
is assigned to ES-ES tunneling, whose excitation process 
involves three charge states. The observed pairs of tun- 
neling lines (c - c', d - d', etc.) always coincide when 
extrapolated to Vp ^ V. More complicated excitations 
are expected for many-electron QDs. Four different peaks 
[e, f, c' and d' in Fig. 2(b)] are observed between the 

= 3 and 4 CB regions. The corresponding transitions 
are indicated in the total energy diagram in Fig. 1(b). 
The fluctuation in the total energy can be significantly 
greater than aVp. Non-equilibrium transport can lead 
to the accumulation of energy in excess of the excitation 
energy. As A" increases, many long-lived ESs with differ- 
ent S can contribute to the transport. The complexity 
of many-body excitations increases with aVp and A" Q . 



measurements, we expect to see related features in con- 
ventional dc excitation measurements. For example [see 
Fig. 4(b)], novel DET involving A" = 1, 2, and 3 is 
allowed in dc measurements when eVb is greater than 
|/7(1, ^) - t/(2,0)| (energy required for normal SET be- 
tween A- - 1 and 2 GSs, arrow i), |/7(2,1) - /7(1,|)| 
(excitation energy to the N = 2 spin-triplet ES from the 
A' = 1 GS, arrow ii), and |/7(3, \) - U{2, 1)| (extra en- 
ergy for novel DET, arrow iii). The necessary conditions 
for the novel DET are satisfied in the red region in the 
eVb-Vg plain of Fig. 4(c). Note that this complex exci- 
tation can only occur if Tspin ^ F"^. Figure 4(a) shows 
the tunneling current spectrum of QD-2 between A" = 
1 and 2 CB regions at 5 = T with no Fc-pulse. In 
addition to the expected current step marked by arrow 
a associated with excitation to the ES(2, 1), we see an 
extra current step marked by arrow b associated with 
the novel DET [|l^. Note that the condition in dc mea- 
surement for novel DET as well as the other tunneling 
processes is equivalent to that in pulse measurements, by 
taking aVp equivalent to eVb. 

In summary, we have discussed novel DET, which can 
lead to considerable charge fluctuation, and ES-ES tun- 
neling, which gives rise to significant fluctuation in the 
total spin. These non-equilibrium transport processes 
cannot be explained by orthodox Coulomb blockade the- 
ory, and arise from the absence of effective spin relaxation 
inside a quantum dot. With miniaturization of semicon- 
ductor devices selection rules for tunneling transitions 
have to be considered when manipulating spin and charge 
in a quantum dot. 
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FIG. 4: (color), (a) Color plot of the dc current / in the Vb - 
Vg plane of QD-2 at B = T. The current increases stepwise 
along the line marked by arrows a and b, which are assigned 
to normal GS-ES tunneling and the novel DET, respectively, 
(b) Total energy U{N, S) for the dc excitation processes, (c) 
Corresponding conditions for the dc excitation transport in 
the eVb - Vg plane. Novel DET appears in the red region. 

Although our findings till now are deduced from pulse 
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